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Warm Extended Dense Gas Lurking At The Heart Of A Cold 

Collapsing Dense Core 
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ABSTRACT 

In order to investigate when and how the birth of a protostellar core oc- 
curs, we made survey observations of four well-studied dense cores in the Taurus 
molecular cloud using CO transitions in submillimeter bands. We report here 
the detection of unexpectedly warm (~ 30 — 70 K), extended (radius of ~ 2400 
AU), dense (a few times 10 5 cm~ 3 ) gas at the heart of one of the dense cores, 
L1521F (MC27), within the cold dynamically collapsing components. We argue 
that the detected warm, extended, dense gas may originate from shock regions 
caused by collisions between the dynamically collapsing components and outflow- 
ing/rotating components within the dense core. We propose a new stage of star 
formation, "warm-in-cold core stage (WICCS)", i.e., the cold collapsing envelope 
encases the warm extended dense gas at the center due to the formation of a 
protostellar core. WICCS would constitutes a missing link in evolution between 
a cold quiescent starless core and a young protostar in class stage that has a 
large-scale bipolar outflow. 

Subject headings: stars: formation — stars: pre-main sequence — ISM: clouds 
- ISM: individual (L1521F, MC27) — submillimeter 
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1. Introduction 

Low mass stars form through gravitational contraction within cold dense condensations, 
called dark cloud dense cores. The dense cores remain roughly isothermal as long as radiative 
cooling is comparable to the heating from compression of the medium. This compression is 
due to the gravitational contraction over ~10 5 years (e.g., Masunaga et al. 1998). When 
the central density reaches 10 10 ~ n cm -3 , compression heating becomes dominant and a 
cond ensed, more or less star-sized object, prot ostellar core, forms at the heart of a dense 



core (IBoss fc Yorkelll995l ; iMasunaga et al.lll9981 ). 



A starless dense core is a cold (~ 10 K) dense (> 10 4 cm~ 3 ) large (~ 20,000 AU) 
condensation of gas and dust (e.g., Bergin & Tafalla 2007). It doesn't harbor any embedded 
infrared point source (hereafter IPS). Some evolved starless cores have dynamically collapsing 
envelopes. On the other hand, class objects (e.g., Andre et al. 2000), which are embedded 
in protostellar dense cores with a temperature of ~ 30 K, a density of > 10 5 cm -3 , and a 
size of < 10,000 AU, are characterized by the existence of IPSs and large-scale molecular 
bipolar outflows driven by protostars. Note that there are dense cores that harbor very 
low luminosity embedded sources within them but do not have large-scale bipolar outflows. 
These objects were used to be categorized as starless cores, however, embedded sources were 
revealed by deep water maser emission survey at the Nobeyama 45-m telescope (e.g., Furuya 
et al. 2001) and deep infrared observations by the Spitzer space telescope (e.g., Dunham et 
al. 2008). They are in a stage more evolved than the starless core stage and younger than 
the class stage. The formation of the protostellar cores must occur in such a transition 
stage. 

In order to investigate when and how the birth of a protostellar core occurs, we observed 
four dense cores in the Taurus molecular cloud using two CO lines in submillimeter bands, 
i.e., J = 6 — 5 (wavelength A 433.544 /mi) and 7 — 6 (A 371.650 /iin) transitions. The energy 
levels at J = 5 and J = 6 correspond to 83 K and 116 K, respectively. The critical densities 
at the temperatures of 100 K traced by the molecule are several times 10 5 cm -3 . With the 
molecular lines that trace both warm and dense gas components, one can study the physical 
properties of the regions closely related to the final processes of protostellar core formation 
when material at the center starts warming up due to the formation of the protostellar core. 



2. Source Selection 



Our surveyed dense cores include L1521E, L1498, L1521F and L1544. These cores 
are selected because (i) they are well studied and known as objects in early stages of star 
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formation, i.e., the starless core stage or the stage of a dense core that harbors a very low 
luminosity IPS and doesn't show a large-scale bipolar outflow, and (ii) they are reasonably 
isolated hence good targets to study internal star formation processes without any external 
disturbance. Among the aforementioned dense cores, L1521E is in the youngest evolutionary 
stage because of its richness in carbon-chain molecu les (e.g., Hirota et al. 2002). All the 



cores except for L1521E show a gas i nfall signature (IHirota et al.ll2002l ; iTafalla et al.lll998 



Onishi et al.lll999l ; ITafalla et al.ll2004j ). L1498 has a relatively low density of 1 4 cm 3 and is 
in a younger evolutionary stage as compared to L1544 and L1521F . L1521F (IMizuno et al. 



r 



1994! ; ICodella et all 119971 : lOnishi et all 1199a IShinnaga et all 120041 ; ICrapsi et all 120041 ) has 
been one of primary targets among dense cor es because of the sh ort dynamical time scale due 
to the high density of 10 6 cm" 3 at the cen ter (lOnishi et al.lll999l ). chemically evolved features 
Jshinnaga et alihooi ICrapsi et all 120041 ) . and the discovery of the IPS, L1521F-IRS, which 
may be a protostar i n the making, within the dense core us ing the Spitzer space telescope 
(ITerebey et al.ll2005l ; iBourke et al.l 12006c iTerebey et al.l 120091 ). L1544 shows a similar nature 
to L1521F, however, no IPS has been detected to date. Based on the observed properties, 
L1521F is the most evolved sample in the set. 

For L1521F, the Spitzer observation imaged a bipolar nebula of scattered l ight extending 



in th e east-west direction associated with L1521F-IRS (IBourke et al.l 12006c ITerebey et al. 



20091 ). which may trace a cavity. However, no large-scale bipolar outflow is observed using 
spectroscopy. A deep centimeter continuum observation confirms that there are no shock- 
ionized inner regions of a bipolar outflow from a ~ 0.1 L Q source (IHarvey et al.l 120021 ) . 
L1521F-IRS is not categorized as a class source owing to lack of a large-scale bipolar 
outflow. It is in an evolutionary category similar to GF 9-2 (e.g., Furuya et al. 2006, 2009) 
that shows very weak water maser emission from an embedded source but does not have a 
large-scale bipolar outflow. 



3. Observations 

The CO observations in submillimeter bands were carried out at the Caltech Submil- 
limeter Observatory (CSO). The CO J = 6 — 5 (frequency v 691.473076GHz, Goldsmith 
et al. 1981) and 7 - 6 [y 806.651801GHz, Schultz et al. 1985) data were taken on 2006 
February 06 — 08 and on 2009 January 22 UT, respectively. The high-altitude dry site, 
sensitive receivers and high efficiency of the telescope combined to permit the observations. 
The Dish Surface Optimization System was used during the observations to compensate for 
the gravitational deformation of the 10.4-m diameter reflector (Leong et al. 2005). The 
telescope's pointing was checked every one - two hours using planets. The pointing accuracy 
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is estimated to be ~ 3" for both observations and the accuracy of the map registration of 
the two CO transitions is within 3". 

Cryogenically cooled SIS receivers operation at 4 K at the CSO that we used for the 
observations produced typical singl e sideband system temperature of ~ 3900 K at 433.5 fim 



and ~ 2300 - 3000 K at 371.7 (im (IKooi et al. 1120001 1 measured with a 50 MHz bandwidth 



spectrometer. Beam chopping method was used. The beam sizes at the two wavelengths are 
about 10"and 9" at 433.5 and 371.7 /zm, respectively. From observations of Mars and Saturn, 
the main beam efficiencies are measured to be ~30 % and ~35 % at 433.5 /im and 371.7 yum 
at elevations of 30 and 48 degrees, respectively. The RMS noise levels for both observations 
are about 0.6 - 1.4 K in T m ^. The velocity resolution of both observations is ~ 0.02 km s _1 . 



4. Results 

Among the observed dense cores, only one, L1521F, shows significant CO line emission at 
both transitions (Figur^U). The rest of the dense cores don't show the CO 6 — 5 line emission 
brighter than 0.3 K in the observed main beam temperatures (T m b). Observed properties of 
CO lines of L1521F are summarized in Table [TJ T m b values of the CO transitions are weak 
(see Table CD), suggesting that the emission may be optically thin. The peak T m b is observed 
at the position of R.A. = 4 h 25 m 34. s 29, Dec. = 26°45'8".7 (B1950.0) for both transitions. 
Note that L1521F-IRS is located at 9" south of the peak position. Considering the thermal 
velocity width of the CO gas at 100 K is only about 0.5 km s _1 , a significant fraction of the 
broad velocity widths originates from non-thermal motions. 

We estimate the temperature of the observed warm dense gas components. Taking 
advantage of the similar beam sizes of the observations of the two transitions, one can derive 
the excitation temperature (T ex ) using the following equation (e.g., Shinnaga et al. 2008): 

38.714 K 

T, 



r 7"\ 2 Tmb CO 6-5 



hi 



(I)' 



CO 7-6 



T ex values estimated at the positions of (R.A. offset, Dec. offset) = (0", 0"), (-10", +10' n 
(-10", 0"), and (-10", -10") are 57 ± 1.7, 40 ± 7, 74 ± 19, and 68 ± 5.5 K, respectively. 
Note that the CO 7 — 6 emission is not detected towards the position of L1521F-IRS, in- 
dicating that the temperature near the IPS is still low (< 30 K). Over the central 30" x 
30" region (corresponds to 4200 AU), the averaged T ex is calculated to be 34.4 ± 0.94 K. 
Th ese temperatures are much higher than the T ex measured with an N2H + transition, 5.0 



K (IShinnaga et all 120041 ) . 
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Figure [2] presents the total integrated intensity map of the CO J = 6 — 5 emission that 
traces the warm extended dense gas (hereafter WEDG), overlaid on the dust continuum 
map at 850 u rn that traces the co ld (~10 K) extended (16,000 AU) dense (> 10 5 cm -3 ) 
condensation (IShinnaga et al.ll2004l ). The observed effective radius of WEDG within the 3 
a contour is 2400 AU (0.012 pc). WEDG sits in the central region of the cold condensation 
and is distributed somewhat asymmetrically, extending to the north and northwest directions 
from L1521F-IRS. The peak T m b of the CO 6 — 5 transition is found at the position of (R.A. 
offset, Dec. offset) = (0", 0"), while the CO 6 — 5 total integrated intensity map peaks at 
the position near L1521F-IRS, i.e., (R.A. offset, Dec. offset) = (0", -10") owing to a large 
line width measured at the position (observed line width Vfwhm of 2.46 ± 0.46 km s _1 ). 



5. Discussion 

We considered four candidate mechanisms for the excitation of WEDG. These are ra- 
diative excitation from L1521F-IRS, infrared pumping from L1521F-IRS, heating by UV 
radiation generated from a possible accretion disk, like the one t hat may explain n arrow 



submillimeter CO line components associated with class objects ( ISpaans et al.lll995l ). and 
heating due to shocks. Only the last one provides a convincing interpretation of the data. 
Radiative excitation from L1521F-IRS is dismissed as L1521F-IRS's luminosity is not high 
enough to pump CO molecules up to such a high J level over > 2400 AU region. In fact, 
CO 7 — 6 emission is not detected near the L1521F-IRS position. Infrared pumping is ruled 
out due to the large vibrational level spacing of CO (~3000 K, Carroll et al. 1981). The 
protostellar core's temperature doesn't become higher than 3000 K as molecular hydrogen 
begins to dissociate at ~ 2000 K. It is implausible that UV radiation from the accretion disk 
heats the observed CO gas as the peak T m b values of both transitions are not found towards 
the position of the IPS. 

O n the other hand, sh ocks would be able to pump the gas up to ~ 100 K over a large 
region iNeufeld et alill995h . In fact, a systematic velocity pattern is not measured in WEDG, 



indicating that WEDG may be excited in shock-heated regions generated by the collision 
between the cold dynamically collapsing components and outflowing/rotating components 
near the center. Considering the presence of L1521F-IRS, it's highly likely that some outflow 
activities have already been initiated. The cavity of the bipolar nebula also indicates the 
existence of a small jet associated with the protostar. Furthermore, L1521F-IRS may have 
a large rotating circumstellar disk with a size of the order of 10 3 AU. The north-south 
elongated feature seen in the 850 fim dust continuum map may indicate the existence of the 
north-south elongated disk that is perpendicular to the cavity axis, i.e., perpendicular to the 
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axis of the bipolar nebula. Based on the measured velocity gradient of 15 km s" 1 pc _1 over 
0.01 pc with a clu mp mass of 0.1 M P) b y using the N 2 H + J = 1 — transition at the center 
of the dense core (jShinnaga et al.l 12004 ). the centrifugal radius of the collapsing medium is 
estimated to be about 10 3 AU. In addition, the asymmetric distribution of WEDG supports 
the view that WEDG comes from the shocks but not from the small jet itself. If the two 
CO lines are in the LTE condition, the lower limit of the mass of WEDG becomes an order 
of 10~ 2 M . 

Figure [3] illustrates a magnified view of the center of the dense core that we propose. 
The cavity with blue color may be in front of the disk, while the cavity with green color 
is behind the disk, along the line of sight. The cavity in front is opening towards the west 
from L1521F-IRS, which may channel part of the outflowing gas, making the shock easier 
to observe on the western side. The north-south extension of WEDG may be due to the 
circumstellar disk elongated along the north-south direction. 



6. Summary and Future Work 

We made survey observations using the two CO transitions in submillimeter bands to 
search for the dark cloud dense cores which harbor newly born protostellar cores and to 
investigate the physical properties of the objects. This study identified a transient stage 
of star formation between the starless core and class stages. As the starless core stage 
progresses, the cold collapsing dense material starts warming up at the center due to the 
formation of a protostellar core and forms WEDG in the central region. We name this new 
stage "warm-in-cold core stage (WICCS)". One should search for WEDG using warm and 
dense gas tracers such as the CO 6 — 5 and 7 — 6 transitions to identify the objects in 
WICCS. This object would constitute a missing link in evolution between a starless core and 
a protostar, yielding an important step toward understanding of the formation mechanism of 
a protostellar core and a protostar. Survey observations of CO transitions in submillimeter 
bands towards dense cores in early evolutionary stages are necessary to add more samples 
in WICCS in order to obtain a complete picture of WICCS. 

This research was performed at the Caltech Submillimeter Observatory, supported by 
NSF grant AST-05-40882 and AST-0838261. This work was also supported by Grant-in- 
Aids from the Ministry of Education, Culture, Sports, Science and Technology of Japan 
(No. 19204020 and No. 20740113). H. S. is grateful to Richard Chamberlin, Brian Force, 
and Hiroshige Yoshida for their support on the observations in 2009 January. 
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Table 1. Observed Properties of the CO Lines 



Transition 


Wavelength 


T mb a 


VFWHM b 






(yum) 


(K) 


(km s -1 ) 


(km s -1 ) 


J = 6-5 


433.544 


1.5 - 5.5 


1.1 - 4.2 


6.77 (3) 


J = 7-6 


371.650 


1.6 - 3.7 


2.5 - 3.5 


6.66 (8) 



a Observed peak main beam temperatures. 
b Observed line widths. 

c Intcnsity-wcightcd mean velocities at the peak intensity position. 
The numbers in parentheses represent 3 a in units of 0.01 km s _1 . 
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Fig. 1. — Profile maps of the CO 6 — 5 (black line) and 7 — 6 (red line) emission of L1521F. 
The central panel (i.e., the origin of the diagram) corresponds to the position at R.A. = 4h 
25m 34.s29, Dec. = 26°45 / 8.7"(B1950.0), where the peak intensities of the two transitions 
are observed. The peak T m b values of the 6 — 5 and 7 — 6 emission are measured to be 5.5 
and 3.7 K with observed line widths Vfwhm of 2.3 and 2.6 km s' 1 by Gaussian fitting, 
respectively. Note that L1521F-IRS is located near the position of (R.A. offset, Dec. offset) 
= (0", -10"). 
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The total integrated intensity map of CO 6 — 5 emission that traces WEDG (white 
contours), overlaid on the grey-scale contou r map of the cold ext ended components traced 
by the dust continuum emission at 850 /xm (IShinnaga et al.l 120041 ) . The contours of WEDG 
are drawn at 5, 7, 9, 11 and 13 a levels (1 a corresponds to 1.1 K km s -1 ). The map origin 
is the same as that of Figure [H The yellow star marks the location of L1521F-IRS. The 
black dots represent the observed positions. The beam sizes of the CO 6 — 5 and 850 /j,m 
continuum maps are 10" and 14", respectively. 
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Fig. 3. — Schematic illustration of the heart of L1521F. The blue arrows indicate the cold 
dynamically infalling medium. The white contour filled with red color represents the lowest 
contour of WEDG shown in Figure El The yellow star marks the location of L1521F-IRS. 
The irregular shaped objects with green and blue colors associated with L1521F-IRS indicate 
the bipolar cavity. The large rotating disk is perpendicular to the cavity axis. The orange 
arrows pointing outwards through the cavity indicate the outflowing gas from L1521F-IRS. 



